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ABSTRACT 

We conducted an exploratory search for quasars at z ~ 6 — 8, using the 
Early Data Release from United Kingdom Infrared Deep Sky survey (UKIDSS) 
cross-matched to panoramic optical imagery. High redshift quasar candidates 
are chosen using multi-color selection in i, z, Y, J, H and K bands. After removal 
of apparent instrumental artifacts, our candidate list consisted of 34 objects. 
We further refined this list with deeper imaging in the optical for ten of our 
candidates. Twenty-five candidates were followed up spectroscopically in the 
near-infrared and in the optical. We confirmed twenty-five of our spectra as very 
low-mass main-sequence stars or brown dwarfs, which were indeed expected as 
the main contaminants of this exploratory search. The lack of quasar detection 
is not surprising: the estimated probability of finding a single z > 6 quasar down 
to the limit of UKIDSS in the 27.3 square degrees of the EDR is < 5%. We find 
that the most important limiting factor in this work is the depth of the available 
optical data. Experience gained in this pilot project can help refine high-redshift 
quasar selection criteria for subsequent UKIDSS data releases. 

Subject headings: surveys — cosmology: observations — quasars: general - 
stars: low-mass, brown dwarfs 
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Introduction 



Quasars remain one of the most powerful probes of the early universe, providing mea- 
surements and constraints for reionization, early structure formation, and early chemical 
enrichment. Absorption spectra of quasars at z > 5.5 provide some of the key measurements 
for the evolution of the inter-gal actic medium (IGM) transmission at the end of th e cosmic 
reionization (see reviews by, e.g. iFan et al.ll2006l ; iDjorgovski et al.ll2006l ; lLoebll2006l . and ref- 
erences therein). This includes direct measurements of the optical depth of IGM absorption, 
and estimates of the radii of quasar Stromgren spheres, which constrain the neutral hydrogen 
fraction in the regimes xhi ~ 10~ 4 — 10~ 2 , and ~ lO^ 1 , respectively. Metallic line absorbers 
seen in their spectra constrain the early chemical evolution of the IGM. Abundance analysis 
of quasar spectra reaching out to z ~ 6 implies a high degree of chemical evolution in their 
hosts, with metallicities up to ~ 10 and Fe/Mg ra tios indicative of a very early chemical 
enrichment by type I SNe (e.g. iDietrich et al.ll2003bl ]a ). They seem to b e powered by black 



holes with masses reaching a few times 10 9 M (e.g., I Vestergaardl 120041 ) . implying an early 
and extremely efficient mass assembly. Luminous quasars may mark some of the densest and 
most biased spots in the early univers e, which may have enhanced rates of galaxy formation 



( iStiavelli et al.ll2005l ; lDjorgovskill2005l . and references therein). 



All of these observations represent non-trivial constraints for theoretical models, which 
become sharper as the quasar redshift increases. The problem in pushing to higher redshifts 
is that the Lya line and the associated flux drop used to color-select quasars move outside 
the CCD wavelength range for z > 6.5 or so. Since the most luminous quasars are extremely 
rare, especially at such high redshifts, a large, moderately deep, panoramic IR sky survey is 
needed. 



The UK Infrared Telescope (UKIRT) Infrared Deep Sky Survey (UKIDSS. lLawrence et al 



20061 ) offers an unprecedented opportunity to discover bright quasars at redshifts beyond the 
z ~ 6.5 barrier which lim its the purely optical surveys such as the Sloan Digital Sky survey 
(SPS S, lYork et al.ll2000l ) and Palomar Quest (PQjj. The luminosity function of (IFan et al. 
20011 ) predicts that about 10 such quasars should be discovered over the 4000 deg 2 , the final 
coverage area planned for the UKIDSS Large Area Survey (LAS). 

With this newly available data, we have embarked upon an exploratory study, using cat- 
alogs from the U KIDSS Early Data Rele ase (EDR) to conduct a preliminary search for high- 
redshift quasars. IVenemans et al.l (120071 ) took advantage of the ov erlap between the UKIDSS 



first data release (DR1) and the SDSS deep survey (Stripe 82; lAdelman-McCarthy et al. 



^ttp : //palquest . org 
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20071 ) which reaches ~ 2 magnit udes deeper than t 



to-near-infrared color selection, IVenemans et all (120071 ) found ULAS J020332.38+001229.2 



re sin gle-scan SDSS images. Using optical- 



az = 5.86 quasar. This illustrates the potential of UKIDSS for finding high redshift quasars. 

The goal of this study is to explore this further, with the available UKIDSS EDR data, 
combined with available optical data from the SDSS and PQ surveys. Our purpose here is to 
understand better the mothodology, the limitations of the available data, and the practical 
issues posed by the selection of high-redshift quasar candidates in this regime. We do not 
try to estimate the quasar luminosity function at these previously unprobed redshifts, since 
the initial available UKIDSS EDR coverage is too small to expect detections of even a few 
z > 6.5 quasars, as we discuss below. 

In §2 we describe our selection technique, including our reliance on SDSS and PQ 
imaging in §2.1 and §2.2, respectively. In §3 we describe our followup observations; we 
present the imaging data in §3.1 and our spectroscopy in §3.2. In §4 we analyze our spectra 
and in §5 we provide a brief discussion of our results. 



2. Selection Technique 



The UKIDSS survey is a near-infrared imaging sky survey which will cover 7500 deg 2 
at its completion. The survey is comprised of five separate surveys with varying depths, 
coverage areas, and filter combinations aimed at specific science goals. The Large Area 
Survey has the ideal design for high redshift quasar searches. Covering ~4000 deg 2 in the Y, 
J, H and K passbands and designed to reach signal-to-noise ratio of 5 at K = 18.1 (Vega 
magnitudes), the survey is primed to detect high redshift (z ~ 7) quasars. 



The UKIDSS Early Data Release (EDR. iDve et al||2006j ) covers 27.3 deg 2 of LAS data. 
The median 5cr point source depths for this data in the four LAS filters are Y = 20.23, 
J = 19.52, H = 18.73 and K = 18.06 (Vega magnitudes). We select our high redshift 
quasar candidates by combining near-infrared data from the LAS in the EDR of UKIDSS 
with optical data from SDSS and PQ an d selecting objects in color-color space based on the 
synthetic colors for quasars presented in iHewett et al.l (120061 ) . 

Using the WFCAM Science Archive (WSA) SQL interfac^] we selected the entire EDR 
catalog by querying for all objects wit hin a set of cones chosen to include the entire LAS 
shown in Figures 9 and 10 of iDye et al.l (120061 ) . We select all objects within a 
90 arcminute radius centered on a series of 17 pointings in the LAS2 and LAS4 regions. 



^ http : //surveys . roe . ac . uk/wsa/edrplus jrelease . html 
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In LAS2 we queried for objects in cones centered on 14 /l 52 m 30 s +06°30', 15 h 29 m +05°40', 
15 h 39 m +06°00', and I5 h 49 m +05°40'. In the LAS4 area, we centered our queries on 13 
equatorial positions (Dec.= 00°00') with R.A .= 12 h 13 m , 12 h 24 m , 12 h 35 m , 12 /l 46 m , 12 /l 57 m , 
^OS" 1 , 13H9 m , 13 h 30 m , 13 fc 41 m , 13 h 52 m , 14 ft 03 m , 14 ft 14 m , 14 ft 25 m . 

We matched the 184245 sources in LAS2 and the 332129 sources in LAS4 to the SDSS 
Data Release 4 PhotoPrimary catalog with a search radius of 2". There were 423214 (151433 
and 271781 matches to LAS2 and LAS4, respectively) matches to SDSS and 93160 (32812 and 
60348 in LAS2 and LAS4, respectively) UKIDSS sources that were unmatched. We required 
that the candidates be classified as 'stellar' or 'probableStar' in the UKIDSS catalog. In 
addition, we require that the optically-detected candidates be classified as stellar in the 
SDSS catalog as well. 

The UKIDSS catalog provides magnitudes for objects measured in a series of circular 
apertures whose diameters range from 1-13". All magnitude s receive a n aper ture correc- 



tion approximating the total magnitude for a point source. iDye et al.l (120061 ) state that 
aperMag3, the 2" aperture, provides the most accurate total flux estimate. We therefore use 
aperMag3 for the four UKIDSS filters when computing the colors of our candidates. Since 
we are searching for point sources, we use the psfMag magnitudes for the SDSS-matched 



ob. 
by 



ects. These a r e the default asinh magnitudes on the AB photometric system as defined 



Lupton et all (119991 ) 



We expect that the objects for which we are searching will be fainter than the SDSS 
detection limit. Therefore, we assign the unmatched UKIDSS sources optical magnitudes of 
i = 22, which corresponds to a signal-to-noise ratio of ~ 1.5. Since undetected objects can 
be many magnitudes below a survey's nominal limit, we assigned the z-band limit to be 22 
magnitudes, which corresponds to < la detection; this is a liberal cut, but one which assures 
that no viable candidates would be missed. We base these signal-to-noise ratio estimates 
on thorough Monte-Carlo simulations of faint artificial sources placed in and extracted from 
SDSS i- and z-band frames. 



2.1. Color Cuts 



Undetected sources in the UKIDSS catalog are listed with a magnitude of -9 . 9999 95E+008. 
To compute colors, we assign such objects the 5a detection limits provide in lDye et al. J2006h : 



Y = 20.23, J = 19.52, H = 18.73, and K = 18.06. The asinh magnitudes utilized by SDSS 
ensures that a magnitude is always computed from the measured flux, regardless how small 
(or negative). We therefore need not assign magnitude limits to these objects. 
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For UKIDSS objects with SDSS matches, we keep only those with J > 17.5 and imposed 
the following color cutg^l : 



(i - z) AB > 2.5 
and 

(z AB - Jvega) < 2.0 and (i AB - Y Vega ) > 3.0 and (Y - J)ve ga < 0.9 

or 

(i - z) AB > 3.5 and (i AB - Y Vega ) > 3.0 and (Y - J) V e ga < 0.9 

or 

(z AB - Jvega) > 3.0 and (J - H) Vega > 0.2 and (J - K) Vega > 0.6 



Assuming the unmatched sources were % and z-band dropouts we set z and i to 22 
magnitudes, apply the lower limits to the UKIDSS magnitudes as for the matched sources 
and keep only objects with J > 18 magnitudes. We then imposed the following color cuts: 



(zab - Jvega) > 3.0 and (J - H) Vega > 0.2 and (J - K) Vega > 0.6 (2) 

or 

(z AB - Jvega) < 2.0 and (i AB - Y Vega ) > 3.0 and (Y - J)ve ga < 0.9 (3) 

The selection criteria in Equations [2] and [3] concentrate on quasar colors in different 
redshift regimes, while avoiding L an d T dwarfs. Based on the quasar tracks and brown 



dwarf colors from lHewett et al.l (120061 ) . Equation [2] selects objects with 5.8 < z < 6.6 using 
the the (z — J) vs. (J — H) and (J — K) color-color spaces, while Equation [3] selects quasars 
with z > 7 using their location in the (i — Y) vs. (Y — J) color-color space. A band of L 
and T dwarfs exists at 2 < (z — J) < 3, which we avoid at the expense of excluding quasars 
in the 6.6 < z < 7 redshift range in our selection. 

After imposing these cuts we had 269 candidates, 2 of which had SDSS detections and 
267 which were undetected by SDSS. We require that the objects with no SDSS counterpart 



3 Using the values directly from their respective catalogs, we present the SDSS-UKIDSS colors as hybrids 
of the AB and Vega magnitude systems. We state the magnitude type explicitly in Equations Q] El and [3J 
but drop the subscripts thereafter. 
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be detected in at least two UKIDSS bands and eliminated single detections, some of which 
were obvious cross-talk artifacts. We also examined 8'x 8' UKIDSS cutouts in all four bands 
for these candidates and eliminated 123 spurious candidates caused by image artifacts. Many 
of these were cross -talk images or other artifacts associated with nearby bright stars (see §7 
of bve et alJbood ). 



2.2. SDSS Images 

Next, we extracted image cutouts for the remaining 146 candidates from SDSS in the r, 
i, and z bands using a Simple Image Access Protocol (SIAP) service available through the 
SDSS Virtual Observatory (VO) toolfl The SIAP allows for the retrieval of SDSS images 
by specifying a position on the sky, an image size, and desired filter. We coadded all images 
containing each candidate which, in some cases, increased the number of passes to three or 
more images, if multiple passes of SDSS are available. In addition, we stacked each source's 
r- and i-band images to assure no flux is detected blueward of z-band dropouts and we 
stacked the r-, i-, and z-band images to look for faint detections in otherwise undetected 
sources. 

We examined these images and eliminated objects affected by image artifacts or proxim- 
ity to very bright stars. In addition, there were several isolated, well-detected sources that 
were not cataloged in the PhotoPrimary catalog, but which were flagged as "secondary" 
SDSS sources. Although secondary objects are defined in the SDSS table schema as re- 
observations of primary objects, these objects did not have a primar y entry in the SPSS 



catalog and were solely secondary (a similar observation was made by lMaddox et al.l 120081 ) 
We removed such objects. 

After these eliminations, we are left with 36 candidates for z > 6.5 quasars, one of these 
has an SDSS match and 35 are unmatched in SDSS. Two pairs of candidates are listed in 
the UKIDSS catalog as separate sources, but are <0.5" apart. We treat each of these as 
a single source, reducing our candidate number to 34. We list their average positions and 
fluxes under a single entry in Table [TJ 

We plot in Figure [1] the optical-near-infrared colors of the candidates which did not 
have SDSS matches. On the left, we plot (i — Y) vs. (Y — J) and on the right we plot 
(z — J) vs. (J — K), using the SDSS limits for the i- and z-band magnitudes. We plot with 
contours all stellar UKIDSS sources with SDSS counterparts, representing the stellar locus. 



4 http: / / skyserver.sdss.org/ vo/dr6siap/SIAP.asmx 
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The candidates that obeyed our color criteria are plotted with small black points, while the 
34 objects that passed our visual inspection are plotted with large filled symbols; the circles 
are objects without counterpart in the SDSS catalog and the square is the candidate with 
a match in the SDSS catalog. We note that although we use two sets of color cuts for the 
two redshift regimes (Equations [2] and [3]) all of our optically unmatched candidates are in 
the higher redshift regime. Only the object that had a match in SDSS is captured by the 
lower-redshift color-selection. 

With the SDSS images in hand, we can individually measure the magnitude limit of our 
candidates in the i- and z-bands, to better position them in the color-color plots. We then 
computed the asinh magnitude of our source by measuring the sky-subtracted counts in a 
4"-diameter aperture and using the zero-point values provided in the SDSS Frames table 
for each frame @. We list these in Table [U When the sky-subtracted flux in the aperture 
is negative then we list the la lower limit magnitude determined from our simulations; 
Zii m = 22.6 and zu m = 21.1. As these measurements show, our estimate of the i- and z- limits 
at 22 magnitudes is generally an underestimate for % and an overestimate for z, though there 
is significant variation. 

To understand the effects of this variation on our object selection we plot in Figure [2] 
the 33 quasar candidates with no SDSS counterpart with square symbols using the now- 
calibrated SDSS magnitudes in the same color spaces as Figure [TJ (i — Y) vs. (Y — J) 
on the left and (z — J) vs. (J — K) on the right. The open squares are spectroscopically 
confirmed low-mass stars or brown dwarfs, while the filled squares are unidentified sources 
(see below). The error bars correspond to the i- and z-band photometric errors from our 
aperture photometry measurements. For comparison, we plot the quasar model tracks with 
three di fferent continuum slo pes from steep to flat, left to right solid lines, respectively, pro - 
vided in iHewett et al.l (120061 ) and a redshifted elliptical galaxy from iMannucci et al.l (120011 ) 
with open triangles. We overplot with aste risks and open circl es L and T dwarfs and M 
dwarfs, respectively, using the tables also in lHewett et al.l ( 120061 ). The filled circle is a very 
cool, T8.5 brown dw arf discovered in the UKIDSS First Data Release (DR1) LAS survey 
( IWarren et al.l 120071 ). We plot this object using our imposed SDSS magnitude limits, as it 
is not detected in SDSS. Both our color criteria (Equations [2] and [3]) reject this object. In 



addition, we overp lot in the right-hand panel with filled triangles five quasars from iFan et al. 
(120001 . 1200 ll . |2003j) for which if -band photometry was available. 



5 This process is described in detail in|http : //www, sdss . org/dr6/algorithms/f luxcal .html 
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2.3. PQ Images 

Palomar- Quest (PQ) is a large, digital, synoptic sky survey, conducted at the Palomar 
Observatory's 48-inch Samuel Oschin telesco pe (P48), using the l arge area QUEST-II cam era, 



a mosaic of 112 CCDs covering 3.6° x 4.6° flBaltav et all 120071 : iDiorgovski et all 120081 ). To 



date, the survey has covered ~ 15, 500 deg 2 with —25° < 5 < +25° in up to four filters, 
either Johnson U, B, R and / or Gunn g, r, i, and z. Typical limiting magnitudes in a single 
pass are ~ 21 magnitudes, and in ~ 6 — 8 coadded passes under good conditions they reach 
the depth of SDSS in the redder bands. The data are taken under almost all reasonable 
conditions, including non-photometric and poor seeing. Image full-width at half maxima for 
stars vary between 2" and < 4". 

Although the survey data have not yet been fully photometrically calibrated, a mere 
detection in an R- or r-band filter is sufficient to eliminate a candidate. We examined co- 
added PQ images of our 34 candidates in Gunn r, i and z, and Johnson R and /. One 
candidate, 132022.09+000211.52, was detected in a Johnson R-band image at 13 fc 20 m 22!08 
+00°02'12'.'53, only 1'.'02 away. The co-added PQ image is made up of eight passes. In 
addition, the Johnson passband is broader and redder than the SDSS Gunn filter reaching 
a deeper flux limit in each scan. The object is visible by eye and a 3" diameter aperture 
measured the source at a 2a r detection. 132022.09+000211.52 is also detected in the PQ 
Johnson /-band (Figure [3]) with a ~ 3.5a/ detection, though a possible weaker detection can 
be seen in the corresponding SDSS image on the right. 



3. Observations 
3.1. Imaging 

We obtained images for ten of the 34 candidates at the 200 inch Hale Telescope at 
Palomar Observatory. Seven fields were observed with the Large Format Camera (LFcfj 
on UT 2006 May 31. The sky was clear and the seeing was ~1". Since the LFC field of 
view is 25.3' x 24.8' several nearby sources were detected on two of the images. All but 
one of the fields were observed with a single 300 second integration. One field, containing 
130248.16—001917.70, was observed in three dithered frames for 120, 120 and 300 seconds, 
which we combined to create a final image. 

Since any detection in the r or % bands is sufficient to eliminate z > 6.5 quasar candidate, 



e- 



http : //www . astro . caltech . edu/palomar/200inch/lf c/lf c_spec . html 
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we used a broad- RI filter with central wavelength of 7670A and width of 2940A to collect as 
many r- and z-band photons as possible in a relatively short integration time. All ten sources 
were detected and therefore eliminated as candidates. Figure S] shows the depth gained by a 
300 second exposure image with the P200 (left) over the combined r+i+z band SDSS image 
(right). 

We reduced the LFC data using the IRAF MSCRED package following standard image 
data reduction procedures, including an astrometric calibration that was necessary in order 
to identify our source in the crowded fields, as these images go ~ 2 magnitudes deeper than 
SDSS. We calibrated th e photometry of our RI observations by extracting catalogs with 



the SExtractor package (IBertin fc Arnoutslll996l ) and matching the positions to the SDSS 



catalog. Choosing ~ 150 well-detected sources in SDSS that were not saturated in our LFC 
images, we derived the photometric solution from the r- and i-band SDSS data. The 5a 
magnitude limit of these images is ~ 24 ±0.2 magnitudes. We measure the magnitude of our 
candidates in a 2" diameter aperture centered on the position of our candidate and report 
it in Table [TJ 

Some of these objects were included in the spectroscopic sample, as indicated in the last 
column of Table [fl All of these were confirmed to be low-mass stars. 



3.2. Spectroscopy 

Spectroscopic followup was conducted in the optical and near-infrared for 27 of our 
candidates. Of these we observed all candidates which had no P200 observations, as well as 
seven of the ten objects with a detection in the i?/-band P200 images. 

We obtai ned spectra of 15 3852.53+053714.43 and 145742.42+064058.96 with the LRIS 



spectrograph (lOke et al.lll995l ) on the Keck I telescope on UT 2006 May 21. We used the 
red camera with the 600 line mm -1 grating blazed at 10000 A corresponding to a wavelength 
coverage of ~ 7900 — 9500 A. This allows for the detection of Lya for z < 6.8. We obtained 
a single 900 second exposure for each object, and although the continuum was detected, 
the strong fringing could not be removed with a single exposure. We show the spectra in 
Figure [5j The spectra are similar in shape and are suggestive of late M to early L dwarf 
spectra. We identify strong stellar absorption features: Na I A A 8183,9195 A doublet, a TiO 
absorption band at 8450A, and Ca II A8542 A. We plot in the bottom panel the OH line 
emission spectrum to identify features from the fringing in the object spectra. The presence 
of a continuum with no Lya detection allows us to rule out the possibility of a z > 6.8 
quasar. 
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Anothe r 29 spectra for 26 objects were obtained with the near-infrared spectrograph 
NIRSPEC ( IMcLean et al.lll998l ) on the Keck II telescope on Mauna Kea during three ob- 
serving runs. We observed during two half-nights on UT 2006 June 11 and June 16, one 
full night on UT 2007 January 24 and two two-hour blocks on UT 2007 February 08 and 
February 09, totaling ~2.5 nights. The weather was clear during all the runs, with the seeing 
< 1". We observed our candidates in the low-resolution mode using a 0.76" slit for all nights 
except on UT 2006 June 11, where the ~ 0'.'5 seeing allowed the use of a / .'58 slit. Most of 
the objects were observed with 1200 second exposures; a few brighter objects were observed 
with 600 and 300 second exposures. 

We used the NIRSPEC-3 ( J-band) filter, which spans 1.143 /xm-1.375 /im, for our initial 
spectroscopy. We chose this wavelength window because our exp ected main contamin ant, 
M, L and T dwarfs, are well-studied in this wavelength regime (e.g. . McLean et al.ll2000l ) and 
contain strong neutral potassium (K I) absorption lines at 1.169,1.177 /im and 1.2432,1.2522 
fim. A sky-subtracted, dithered pair of well-exposed spectra immediately revealed the K I 
absorption features, making for a very efficient observing program. On all nights, except 
UT 2006 June 11, we observed at least one A0V star at an airmass representative of the 
airmasses of our target observations. We were unable to obtain a calibration spectrum on 
UT 2006 June 11, and therefore calibrated these data with the A0V star from June 16. 

The data were reduced using the IRAF WMKONSPEC package designed for reducing 
low-resolution NIRSPEC datg@. The two-dimensional spectra produced by NIRSPEC are 
curved and distorted with respect to the detector rows. This software removes the spatial 
distortions allowing for proper extraction, sky subtraction, and wavelength calibration. We 
corrected for telluric abs orption in th e se spe ctra using the spectrum of the A0V star using 
the technique outlined in IVacca et al.l (120031 ). 



4. Results 

Given the relatively small area coverage of the EDR, it is not surprising that we do 
not find any high-redshift quasar detections. In this exploratory search, we have made 
deliberately inclusive and liberal color selection of candidates so as to assure that no real 
quasars would be missed. It is expected that most of the photometry contaminants in this 
parameter space would be stars at the bottom of the main sequence or brown dwarfs. This 
is indeed what we found in our spectroscopic followup. 



http : / / www2 . keck . hawaii . edu/inst/nirspec/wmkonspec/ index . html 
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Figure [6] shows five of our better-detected J-band NIRSPEC spectra chosen from the 
low-mass stars or brown dwarfs found in our survey. The vertical dotted lines indicate the 
position of K I absorption lines at 1.1690, 1.770 jum and 1.2432,1.2522 /im. To help identify 
spectral features that may result from poor sky subtraction or calibration, we plot in the 
bottom panel of each column the mean OH line emission spectrum and the mean telluric 
absorption spectrum for our data. This is especially helpful for identifying spectral features 
in poorly calibrated spectra. 



In an effort to classify these stars, following McLean et al.l (120001 ) and McLean et al. 



(120031 ). we measured the equivalent widths (EWs) of the four K I absorption lines as well 
as the strength of the H 2 absorption at ~ 1.34 /mi. To determine t he line EWs, we fit a 
multiple component Gaussian plus linear continuum to each doublet. iMcLean et al.l (120001 ) 
use the 1.33 /mi/ 1.27 /im ratio to measure the strength of the H 2 absorption, because of 
the strong sky emission at 1.27 /im, we used the ratio of the mean fluxes in the 1.335-1.345 
jum and 1.295-1.305 /im ranges in our measurement. Table [2] lists the results, as well as the 
EWs for four comparison objects from IMcLean et al.l (120031 ) . The EWs of these stars are 
significantly smaller than those of L and T dwarfs whose EWs rise from 6 A for early-L to 



> 10 A in late-L and early-T dwarfs ( McLean et al.ll2000l ; IMcLean et al.ll2003l ; ICushing et al. 
2005 ). M dwarfs have K I EWs of ~ 3.5 — 7.5A for spectral types ranging from M6 to M9; 
the median EWs in our sample is ~ 3.5A. In addition, H2O absorption in L and T dwarfs 
is very strong, with 1.33 /im/1.27 /im < 0.55, while our flux ratios are mostly > 0.7. We 
note, however, that this diagnostic is especially unreliable in our spectra because of the poor 
calibration of some sources leading to an artificial rise in the continuum beyond 1.33 /im. 

We were unable to classify three of our objects from their J-band NIRSPEC spectra. The 
K I absorption lines and the expected water absorption beyond 1.33 /tm were not detected 
in the spectra of 131002.98+000751.88 and 133048.43-002316.79 observed on UT 2006 June 
11, where no calibration star was available. We re-observed these objects on UT 2007 Feb 
8, calibrating on an A0V star observed less than an hour after the target observations at 
an airmass difference < 0.01. 154231.67+053205.29 was observed at the end of the night as 
the sky was brightening. Our spectrum has several saturated sky lines which are obvious in 
the spectrum. Nevertheless the expected K I absorption lines at 1.2432,1.2522 /im are in a 
region of weak sky emission and maximal transmission and are absent from our spectrum. 

It is very unlikely that these objects are high-redshift quasars. The only redshift range 
beyond z = 5.7 (where Lya moves out of the i-band) that is absent of emission line features 
in the observed J-band is z ~ 6.1 — 6.4. The colors of these objects, however, do not fit this 
redshift range in either color-color space shown in Figure [2] and are inconsistent between the 
two panels for any redshift. 
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It is possible that these are compact elliptical galaxies at z ~ 1.5. The seeing in 
UKIDSS images can be as large as I'M, corresponding to ~ 11.5 kpc at z — 1.5. In addition, 
133048.43-002316.79 is classified as 'probableStar' in the UKIDSS LAS catalog. These ob- 
jects would require higher resolution imaging to rule them out as galaxies on a morphological 
basis. 

Figure [2] shows the expected colors of an elliptical galaxy and agrees quite well with 
our candidates in the (i — Y) vs. (Y — J) color-space. However, our candidates are too 
blue in (J — K) and/or too red in (z — J) to fit the expected colors of a z ~ 1.5 elliptical 
galaxy. It is even less l ikely t hat these objects are young, dusty starbursts at similar redshifts. 



Pozzetti fc Mannuccil (120001 ) showed that these objects have even redder (J — K) colors than 
old ellipticals. A broader-wavelength ranged spectrum will likely reveal the nature of these 
objects, though they are highly unlikely to be quasars. 



Discussion 



The space density of z ~ 6 quasars with M1450 < —26.7 is 6.4 ± 2.4 x 10 10 Mpc 



assuming H = 65 km s 1 , Q\ = 0.65 and f) p = 0.35 (IFan et al.l 120041 ) . The quasar density 



evolves with redshift as log[p(,z:)] = — 0.51,2 (IFan et al.ll200ll ). at least out to z ~ 6, and we 



can tentatively extrapolate this trend into the as-yet-unprobed z > 6.5 regime. Our selection 
criteria are sensitive to quasars with 6.0 < z < 6.6 and z > 7.0 (see Figure [2]). From this, 
we compute the expected space density of quasars with M U50 < —26.7 and 7 < z < 8 to be 
~ 6 x 10~ n Mpc -3 . The comoving volume of our survey's 27.3 deg 2 in that redshift range 
is ~ 2 x 10 s Mpc 3 . Although we found no candidates satisfying the low redshift criteria, the 
space density of quasars with 6.0 < z < 6.6 is ~ 1.4 x 10~ 10 Mpc~ 3 over a comoving volume 
of ~ 1.4 x 10 8 Mpc 3 . Therefore, the probability of finding a quasar in the lower-redshift 
regime is ~ 0.02 and in the higher-redshift regime is ~ 0.01, given all of the assumptions 
stated above. The total probability for finding a quasar in our survey of the UKIDSS EDR 
is 0.03. Thus, we did not expect to find any such quasars in this exploratory study, but a 
detection would have been both interesting and significant. 



It is possible that this probability is higher for UKIDSS-selected quasars, since lVenemans et al 



( 120071 ) found, using color selection criteria very similar to those employed here, ULAS 
J020332.38+001229.2, a z = 5.86 quasar in a subset of the UKIDSS first data release (DR1) 
covering ~ 4 times larger area and overlapping the SDSS Stripe 82 which reaches ~ 2 mag- 
nitudes deeper than the single-scan SDSS images. Its absolute magnitude of M1450 = —26.2 
is half a magnitude fainter than the limit for SDSS-selected high-redshift quasars. This 
demonstrates that searches for high-redshfit quasars using UKIDSS data and the methodol- 
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ogy described here are certainly plausible, even at this early stage. 



Jiang et al.l (120071 ) found five z ~ 6 quasars in the SDSS stripe 82 region, including 



ULAS J020332. 38+001229. 2, and determine the space density for quasars brighter than 
M1450 < —26.7 to be 5.0 ± 2.1 x 10~ 9 Mpc -3 , nearly an order of magnitude higher than at 
M1450 = —26.7. At these fainter flux limits, the probabilities of finding such quasars rises 
to 0.15 in the 6.0 < z < 6.6 redshift range and 0.09 in the 7.0 < z < 8.0 redshift range. 
This raises the probability to ~ 25%. However, despite the increase in sensitivity gained 
from selecting for quasars in the near-infrared, our small coverage area still makes finding a 
z > 6.5 quasar unlikely. 

Existing optical datasets that go deeper than SDSS may be useful for reducing stellar 
contamination, as well as probing deeper i nto the high-redshift quasar luminosity fun ction, 



where the density of quasars is higher (e.g. McGreer et al. 2006 ; Mahabal et al.lEood ). 



The m ost important con c lusion of this investigation, as already noted by iHewett et al. 



(120061 ) and IVenemans et al.l (120071 ). is that the searches for high-redshift quasars using 
UKIDSS data will be limited mainly by the depth of the available optical surveys, here 
SDSS and PQ. Since we are selecting objects that are undetected in SDSS, the key limiting 
factor for an expanded highest-redshift quasar search is the depth of the optical imaging. 
We have shown, in Figure HI that a 5-10 minute exposure with the P200 telescope can reach 
1.5-2 magnitudes deeper than SDSS and will eliminate any low-mass star or brown dwarf 
candidates. This is especially important for z > 7 quasar candidates, where no optical light 
is expected and a simple detection in a broad-RI filter image will eliminate the candidate 
from the sample. It is possible that in order to fully exploit the potential of UKIDSS for 
discovery of high-redshift quasars, the existing panoramic optical surveys will need to be 
supplemented by additional wide-field imaging in the far red bands. 

An additional area of possible improvement is the computation of updated sets of quasar 
redshift tracks in these color spaces, which would include both measureme nts of the optica l 
depth to Lya and Ly/3 photons as a function of r edshift beyond z ~ 5 (IFan et al.l 120061 ). 



and a spread of quasar UV-continuum shapes (e.g., iRichards et al.ll2003l ). This can be used 
to further optimize the color separation of quasar candidate from very red stars and brown 
dwarfs and provide a better understanding of the high-redshift quasar selection function. We 
are currently exploring these issues, and will report on their application in a future paper, 
although we do not expect a major difference in the resulting color-criteria. 



We thank Dr. K. Cruz for lending her expertise in low-mass stars to this paper. We 
are grateful to the staff of W. M. Keck observatory for their assistance during our remote 
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Table 1. Color-Selected Candidates 



R.A. 


Dec. 




i 


z 


Y 


J 


H 


K 


RP 


Notes 


(J2000) 


(J2000) 




(AB mag) 


(AB mag) 


(Vega mag) 


(Vega mag) 


(Vega mag) 


(Vega mag) 


(AB mag) 




i n p- o jo ci n 

12 59 49.89 


i o o o 1 1 o 

+00 01 13 


.94 


nn p* o i o a n 

22.50+0.42 


oi no i o p* o 

21.23+0.59 


19.59 


18 


.83 


18.16 


17.49 




b d 


1 O f\ f\ O 1 T J 

13 00 31.74 


i on 1 1 oo 

+00 11 08 


.63 


nn f 1 i o a o 

22.61+0.40 


oo on i o jo 

20.92+0.48 


19.74 


18 


.48 


17.74 


17.06 




b 


1 o m o o pt j 

13 01 39.74 


OO O ^7 i o 

—00 07 13 


.70 


no m i o f o 

23.47+0.53 


OO ITP' I O 1 o 

20.75+0. 18 


20.41 


19 


.52 


18.73 


18.06 






1 O f\ C\ JO 1 

13 02 48.16 


OO 1 O 1 PT 

—00 19 17 


,70 


>22.6 


OO OO 1 O O A 

20.83+0.34 


19.99 


18 


.88 


18.01 


17.29 


nn a a i o oo 

22.44+0.09 


b 


1 O AO OO r* J 

13 03 03.54 


i OO 1 /"> C\ 

+00 16 27 


.66 


no on i o oo 

23.32+0.80 


oo i n i o npr 

20.42+0.27 


20.21 


19 


.00 


18.73 


18.06 






13 03 48.63 


i oo 1 1 4 r 

+00 ll 45 


.02 


>22.6 


ni opt i o r* o 

21.07+0.58 


19.76 


19 


.00 


18.09 


17.73 


nn r> f i o i o 

22.65+0.13 


b 


13 06 33.64 


i O O O O A A 

+00 08 44 


.84 


22.55+0.40 


no 1 o i o -~\ • 

20.18+0.27 


19.66 


18 


.80 


18.39 


18.06 






13 09 55.32 


■ o o n p* n a 

+00 25 24 


.21 


no pro ii c~\ a 

23.78+1.04 


n 1 n i /~\ pto 

21. 86+0. 70 


19.48 


18 


.57 


18.15 


17.49 




b 


13 10 02.98 


+00 07 51 


.88 


n n o o i o o o 

22.30+0.30 


n i no i o /~* a 

21.23+0.64 


20.33 


18 


.97 


18.31 


17.91 




b 


I o no o n a a 

13 20 22.09 


i OO O O 11 

+00 02 ll 


.52 


n 1 o o i o n a 

21. 92+0. 24 


oi no i o r" o 

21.29+0.59 


19.52 


18 


.74 


18.22 


17.62 




b 


1 o no or" o o 

13 28 35.98 


oo on i o 

—00 22 13. 


.99 


nn r* f i o j t 

22.65+0.47 


oo opt i o no 

20.37+0.23 


19.58 


18 


.75 


18.18 


17.52 


nn a a i o 1 1 

22.44+0.11 




■ . t no p" J p" 

13 28 54.53 


oo no f 1 

—00 23 51. 


.32 


>22.6 


n 1 a n i o /■» p* 

21. 42+0. 65 


19.78 


18 


.80 


17.94 


17.26 


nn pt r> i o o 1 

22. 76+0. 31 


b 


1 o on jo jo 

13 30 48.43 


oo no i 

—00 23 16 


.79 


ni oo i o of 

21. 92+0. 25 


oi on i o p* n 

21.32+0.52 


19.96 


18 


.59 


17.47 


16.69 




b 


■ o on 1 r\ o r\ 

13 32 19.39 


i oo no a p* 

+00 20 45 


.43 


nn oo i o opt 

22.39+0.37 


1 o o /~* i o no 

19. 96+0. 20 


19.45 


18 


.72 


17.86 


17.50 




b 


1 o j n p- 1 oo 

13 42 51.33 


i oo no a n 

+00 20 42 


.84 


n 1 o j i o nn 

21. 84+0. 22 


o i on i 1 11 

21.92+1.11 


19.52 


18 


.87 


18.20 


17.82 




b 


10 A A A A -1 A 

13 44 44.14 


oo no r* n 

—00 23 59. 


.36 


no v o i o o j 

23.58+0.94 


ni pt i o f o 

21.67+0.68 


19.98 


18 


.68 


18.06 


17.52 




b 


1 n p* n o i p- n 

13 52 34.52 


i oo i /"i r" o 

+00 16 50 


.43 


nn 1 j i o o o 

22.14+0.30 


no n j i o n p* 

20.24+0.25 


19.36 


18 


.04 


17.01 


16.21 




b 


14 13 32.07 


-00 16 14. 


.84 


24.04+0.98 


20.79+0.34 


20.24 


18 


.99 


18.10 


17.08 




b 


14 18 59.21 


+00 01 28 


.77 


22.33+0.34 


20.74+0.38 


19.24 


18 


.25 


17.47 


17.09 




b 


14 20 08.94 


+00 14 14 


.71 


21.86+0.25 


20.89+0.45 


20.28 


18 


.90 


18.21 


18.06 






14 57 42.42 


+06 40 58 


.96 


21.70+0.21 


20.76+0.33 


19.67 


18 


.70 


18.13 


17.69 




c 


14 57 45.72 


+06 36 29 


.70 


23.7lil.l0 


21.79+0.67 


19.97 


18 


.97 


18.34 


17.32 


22.44+0.10 f 


b 


15 26 51.74 


+05 46 56 


.25 


23.79+1.16 


21.80+0.76 


19.77 


18 


.69 


18.21 


17.61 




b 


15 27 32.25 


+05 20 40 


.81 


22.67+0.47 


21.36+0.76 


19.88 


18 


.71 


17.98 


17.35 


23.15+0.16 


b 


15 27 57.36 


+05 29 14 


.41 


21.94+0.29 


20.81+0.48 


19.52 


18 


.79 


18.21 


17.70 


21.77+0.06 


d 


15 28 03.56 


+05 24 59 


.52 


21.98+0.29 


20.35+0.33 


19.72 


18 


.67 


18.08 


17.42 


22.33+0.10 


b 


15 34 22.57 


+05 45 00 


.62 


>22.6 


>21.1 


19.70 


18 


.43 


17.79 


17.26 


22.33+0.07 


b 



Table 1 — Continued 



R.A. 


Dec. 


i 


z 


Y 


J 


H 


K 


RI a 


Notes 


(J2000) 


(J2000) 


(AB mag) 


(AB mag) 


(Vega mag) 


(Vega mag) 


(Vega mag) 


(Vega mag) 


(AB mag) 




15 37 39.52 


+05 40 10.57 


23.70+1.08 


20.87+0.37 


20.06 


18.72 


18.16 


17.44 


22.77+0.11 




15 38 52.53 


+05 37 14.43 


22.71+0.64 


21.96+0.87 


19.69 


18.77 


18.35 


17.90 




b c 


15 42 31.67 


+05 32 05.29 


21.69+0.22 


21.00+0.40 


19.72 


18.76 


18.00 


17.17 




b 


15 47 00.58 


+05 47 38.03 


23.41+0.82 


20.49+0.27 


19.45 


18.56 


17.97 


17.70 




b 


15 47 21.08 


+05 25 49.45 


22.58+0.52 


21.30+0.55 


19.41 


18.62 


17.85 


17.35 




b 


15 47 27.87 


+05 24 51.90 


23.14+0.79 


20.80+0.36 


19.78 


18.58 


17.66 


16.76 




b 


15 51 27.70 


+05 28 39.07 


22.65+0.55 


20.74+0.36 


19.89 


18.90 


18.67 


17.96 




b 



a Objects in this column were imaged at the 200 inch Hale Telescope at the Palomar Observatory with the Large Format Camera (LFC) 
using a Broad- RI filter, centered at 7670A(AA = 2940A). Integration times were 300 seconds, except 130248.16— 001917.7 which was exposed 
for 540 seconds. 

b This object has a NIRSPEC spectrum. 

c This object has an LRIS spectrum. 

d This object appeared as a pair of candidates separated by < 0'.'5 in the UKIDSS catalog. We treat it as a single source and list its average 
position and flux. 

c This object was selected as a candidates from the list of UKIDSS sources with SDSS counterparts. The i and z magnitudes are the psf Mag 
and psf MagErr entries in the SDSS PhotoPrimary catalog. 

f The magnitude derived for this source is highly uncertain due to a grazing cosmic ray at the location of the source. 
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Table 2. K I Line Equivalent Widths and the Strength of the H 2 Absorption at 1.34 /im 



Wavelength 



Object 


1.169 /Jin 


1.177 iim 


1.244 /im 


1.253 /j,m 


1.34 /im/1.29 




(A) 


(A) 


(A) 


(A) 




1259+0001 


7.6+2.2 


0.8+0.9 


2.4+0.4 


1.8+0.5 


0.90 


1300+0011 


9.4+1.4 


1.1+0.5 


3.8+0.3 


3.3+0.3 


0.61 


1302-0019 


1.7+0.5 


3.9+0.6 


0.0+0.4 


3.7+0.5 


0.71 


1303+0011 


3.4+0.4 


3.7+0.4 


4.4+0.5 


5.7+0.4 


0.88 


1309+0025 


3.2+0.5 


0.8+0.2 


1.8+0.2 


2.6+0.3 


0.76 


1320+0002 


2.9+0.4 


1.6+0.3 


4.6+0.3 


1.2+0.2 


0.88 


1328-0023 


3.8+0.6 


6.7+0.8 


4.3+0.4 


3.0+0.4 


0.36 


1332+0020 


4.1+0.7 


2.8+0.5 


8.6+0.6 


1.9+0.1 


0.80 


1342+0020 


1.7+0.2 


1.8+0.3 


1.2+0.2 


1.9+0.1 


1.23 


1344-0023 


3.6+0.6 


1.6+0.3 


1.7+0.1 


1.7+0.2 


0.69 


1352+0016 


3.4+0.5 


2.2+0.4 


4.3+0.3 


1.2+0.1 


0.70 


1413-0016 


2.8+0.5 


3.1+0.5 


3.2+0.3 


3.3+0.3 


0.73 


1418+0001 


2.5+0.5 


2.1+0.3 


3.5+0.3 


0.6+0.1 


0.74 


1457+0636 


2.6+0.4 


3.0+0.4 


4.5+0.5 


3.1+0.3 


0.70 


1526+0546 


1.7+0.2 


2.4+0.3 


2.9+0.3 


2.8+0.3 


0.70 


1527+0520 


0.8+0.5 


3.3+0.5 


6.0+0.4 


1.0+0.2 


0.81 


1528+0525 


3.1+0.4 


2.8+0.4 


3.8+0.3 


1.2+0.1 


0.74 


1534+0545 


4.9+0.7 


7.0+1.0 


5.0+0.3 


5.4+0.7 


0.83 


1547+0547 


3.1+0.6 


1.6+0.4 


3.7+0.4 


3.7+0.3 


0.70 


1547+0525 


2.9+0.3 


2.0+0.2 


1.8+0.1 


2.6+0.2 


1.35 


1547+0524 


1.0+0.5 


2.9+0.6 


4.3+0.3 


4.7+0.6 


0.67 


1551+0528 


1.5+0.3 


0.6+0.3 


2.1+0.2 


2.1+0.2 


0.81 


M6 - Gl 283B 


3.2+0.2 


4.9+0.2 


3.5+0.2 


3.2+0.2 




M7 - VB 8 


4.4+0.2 


6.4+0.2 


4.6+0.2 


4.6+0.2 




L2 - Kelu-1 


3.7+0.9 


9.0+0.7 


6.9+0.4 


6.8+0.4 


0.55 


L8 - Gl 337C 


5.8+0.3 


8.0+0.3 


3.2+0.3 


4.8+0.3 
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Note. — The last four objects are low-mass stars from iMcLean et al.l (120031 ) whose 
equivalent wid ths are provided for comparison. The 1.34 //m/1.29 /j,m ratio listed for 
Kelu-1 is from IMcLean et al.l ((2000) 
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-1.0 -0.5 0.0 0.5 1.0 1.5 2.0 -1.0 -0.5 0.0 0.5 1.0 1.5 2.0 

(Y-J) (J-K) 

v y Vega v y Vega 

Fig. 1. — Two-color diagrams for our quasar candidates. On the left, we plot i — Y vs. Y — J 
and on the right we plot z — J vs. J — K. The contours represent the density of stellar 
UKIDSS sources with SDSS counterparts (ie., the stellar locus). The small black points 
represent the 267 candidates that obeyed our color selection criteria, including candidates 
that were subsequently eliminated as various artifacts (as described in the text), while the 
filled circles represent the 33 objects with no match in the SDSS catalog that passed our 
visual inspection. The filled square is the only object with a match in the SDSS catalog 
that passed the color criteria in Equation [1] and our visual inspection. The black solid 
line in the right-hand panel corresponds to the color selection criteria [z — J) > 3.0 and 
(J — K) > 0.6 (Equation [2]). The dash-dot line in both panels corresponds to the selection 
criteria (z - J) < 2.0 and (i — Y) > 3.0 and (Y — J < 0.9) (Equation [3]). 
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Fig. 2. — Two-color diagrams for the quasar candidates (squares) using the measured asinh 
magnitudes on the SDSS system. The open squares are spectroscopically confirmed low 
mass stars, while the filled squares are sources with unidentified spectra. The error bars 
correspond to the i- and z-band photometric errors from our aperture photo metry mea- 
surements. The solid lines show quasar model tracks (from iHewett et al.l 120061 ) with three 
different continuum slopes from steep to flat, left to right, r espectively. We plot w ith con- 
nected triangle s the colors of redsh ifted Elliptical galaxy from iMannucci et al. (colors 
obtained from IHewett et al.ll2006T). The as terisks and open circles are L and T dwarfs and 
M dwarfs, respectively ( Hewett et ali 2006 ). The open triangle is the z = 5.9 quasar found 
in th e UKIDSS survey by IVenemans et al.l (120071 ) and the filled triangles are five quasars 
from iFan et al.l J2OO0I . I2OO1L I2OO3I ) for which if - band photome t ry wa s available. The filled 
circle is the very cool T8.5 dwarf discovered by I Warren et al.l (120071 ) in the UKIDSS LAS 
DR1, which is rejected by our color criteria. 
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Fig. 3.— 2'x2' images of 132022.09+000211.52. Left - PQ Johnson /-band image which 
shows a ~ 3.5<7 detection. Right - SDSS z-band image showing weak detection at the location 
of the source. 




Fig. 4.— 2'x2' images of 153739.52+054010.57. The P200 image (left) shows the depth 
gained from a 300 second exposure over the combined SDSS r+i+z image (right). 
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Fig. 5.— LRIS spectra of 153852.53+053714.43 and 145742.42+064058.96 showing stellar 
absorption features, Na I AA 8183,9195A doublet, a TiO absorption band at 8450A, and 
Ca II A8542A. We plot in the bottom panel the mean OH line emission spectrum for easy 
identification of sky lines in our spectra. 
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Fig. 6. — A sample of our better-detected J-band NIRSPEC spectra of stars found in our 
survey. The vertical dotted lines indicate the position of K I absorption lines at 1.1690,1.770 
/im and 1.2432,1.2522 /iin. We plot in the bottom panel of each column the mean OH line 
emission spectrum and the mean telluric absorption spectrum for our data. 



